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Assuming that a quantum field theory with a ^-vacuum term in the action shows non- 
trivial ^-dependence and provided that some reasonable properties of the probability distri- 
bution function of the order parameter hold, we argue that the theory either breaks spon- 
taneously CP at 8 — 7r or shows a singular behavior at some critical C between and n. 
This result, which applies to any model with a pure imaginary contribution to the euclidean 
action consisting in a quantized charge coupled to a phase, as QCD, is illustrated with two 
simple examples; one of them intimately related to Witten's result on SU(N) in the large 
' N limit. 

: 

O ■ 

■ §1. Introduction 

o : 

£^ | Quantum Field Theories with a topological term in the action are a subject of 

interest in high energy particle physics and in solid state physics since a long time. In 
particle physics, these models describe particle interactions with a CP violating term. 
The extremely small experimental bound for the CP violating effects in QCD (strong 
CP problem) is still waiting for a convincing theoretical explanation^ In solid state 
physics, chains of half-integer quantum spins with antiferromagnetic interactions are 
related to the two-dimensional 0(3) nonlinear sigma model with topological term at 
^ ■ 9 = 7T. It has been argued that this model presents a second order phase transition 

at 9 = t:, keeping its ground state CP symmetric (Haldane conjecture).^ 

Non-perturbative studies of field theories with a vacuum term are enormously 
delayed because of two technical reasons: first, there are difficulties in finding a con- 
sistent definition of topological charge for QCD in discrete space-time, preventing the 
use of the most powerful non-perturbative method, the lattice regularization; second, 
even in simpler QCD-like models, as CP N or sigma models, for which a consistent 
definition of topological charge is known, the complex nature of the euclidean action 
forbids the application of all standard Monte-Carlo algorithms to perform numerical 
simulations (see for example Refs. |3J) and[lj)). We need therefore new ideas which, 
beside our present technical abilities, allow a progress in this field. 

The #-term breaks explicitly a Zi symmetry (CP), except at 9 = and 9 = 
ir. Most of the models the 9 dependence of which is known break this symmetry 
spontaneously at 9 = tt (for an early study see Ref. EJ)). An intriguing question in 
such a case is the realization of the 2tt periodicity in #P For instance, in SU(N) 
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gauge theory at large N the 9 dependence appears through the combination 9/N. 
The periodicity is realized by conjecturing that the vacuum energy density is a multi- 
branched function of the form 



E{9) 



N 1 



min H 

k 



9 + 2nk\ 



N 



(1-1) 



This behavior occurs also in some two dimensional systems JZUHIJ3 an( j leads to 
spontaneous breaking of CP at 9 = n in a natural way. Recently, Witten has shown 
that in SU(N) gauge theories at large N Eq. (fT~Tj) holds with H(9) = C9 2 , where C 
is a positive constant independent of TVP^J 

We want to show in this paper that the above behavior is not specific of large iV 
gauge theories but rather general. More precisely, we will argue that any model with a 
quantized charge which is an order parameter of a given symmetry and which appears 
as an imaginary contribution to the euclidean action, either breaks the symmetry at 
9 = it or shows some singularity for 9-values between and it, provided that there is 
a nontrivial 9 -dependence. 

§2. General description 

In the specific case of QCD, the partition function with a topological term in 
the action reads as follows: 



[dA a u ] [#] [#] exp 



d 4 x [C(x) - i9X(x)j 



(2-1) 



16-7T 



7 ^fiupcr Tr F^ v Fp a {x) 



where C(x) is the standard QCD Lagrangian and X{x) 
is the euclidean local density of topological charge, the normalization of which has 
been chosen in such a way that the topological charge, J d 4 x X(x) be an integer. We 
will assume in what follows that the regularized theory preserves the quantization 
of the topological charge. 

To start with the proof, let us write the partition function Zy{9) in a finite 
space-time volume, V, as a sum over all topological sectors, labeled by the integer n 
that gives the topological charge, of the partition functions of each sector, weighted 
with the proper topological phase: 



z v(P) = ^2 9v(n)e 



id ii 



(2-2) 



where 



9v(n) 



exp 



d x C(x) 



(2-3) 



is the standard partition function computed over the gauge sector with topological 
charge equal to n. The function gv( n ) is, up to the normalization factor ^] n gy(n), 
the probability pv{n) of the topological sector n at 9 = 0. If we define the mean 
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topological charge density as x n = n/V, we can write the previous partition function 
as 

Z v {9) = Y, h v (x n )e ieVx -, (2-4) 

where hy{x n ) = gv( n ) and the step Ax n in the topological charge density is 1/V. 

We will try now to write the partition function Zy{6) as a sum of "continuous" 
pseudo-partition functions, where the word continuous stands to indicate that the 
density of topological charge x entering the definition of this new pseudo-partition 
function will be a continuous variable. To this end, let the new hy(x) be any con- 
tinuous smooth interpolation of hy{x n ) and let us define a new function of in the 
following way: 

Z cy (6) = f dxhy{x)e i6Vx . (2-5) 

Summing up the pseudo-partition functions Z c y(6 + 2Trm) for all integers m and 
using the following representation of the periodic delta function (Poisson summation 
formula): 



m m 

we get the following identity: 

Z v (9) = V Y z c,v(0 + 2vrm) , (2-7) 

m 

which relates our QCD partition function Zy{9) to the pseudo-partition functions 
Z c ,v{Q + 27rm). Of course different interpolations of hy{x n ) will give different equiv- 
alent realizations of Eq. Q2-7|) . 

In deriving the above result we have assumed that the sum in m and the integral 
in x can be commuted. This is a working hypothesis on the properties of the hy(x n ) 
function and we will use it to derive some physical consequences. As will be evident in 
the next section, the validity of this assumption is confirmed in some examples where 
analytical calculations are possible. Notice also that the sum in ()2-7|) has nothing to 
do with the sum over different topological sectors. Each sector labeled by an integer 
m in (|2-7j) can get contributions from all topological sectors. If Z c y{0 + 2irm) 
is positive for any integer m, the physical interpretation of these integer numbers 
is very simple. Indeed in such a case and in what concerns the infinite volume 
limit, the summation in (|2-7j) can be replaced by the maximum in m of Z c y{9 + 
2-7rm) (saddle point solution) and therefore the integers m will define metastable 
vacua which will decay to the true vacuum, the energy density of which is given by 
min m {— 1/V log Z c y(6 + 27rm)} (see the first example discussed below). 

The mean value of the density of topological charge can be written as 

(x) ~ V (x) Z c,v(° + ^rn) 
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with 




/ dxxh v (x)e^ 0+27rr ^ V: 
f dxh v (x)e i ( 9+2nm ^ Vx 



(2-9) 



Since x in Eqs. (|2-5|) and (|2-9|) is a continuous variable, the pseudo-partition function 
Z c y{8) wm n °t be a periodic function of 8 and (x) c ,m needs not to vanish at 8 = 
tt. The periodicity in 8 of Zy{0) is however guaranteed because we must sum up 
Z c ,v{8 + 2-Km) for all integer values of m. In the same way, CP is recovered as a 
symmetry of the action at 8 = tt, because the m = contribution to (x) in (|2-8|) 
is, at = tt, compensated by the m = —1 contribution, the m = 1 contribution 
compensates with m = —2 and so on. Therefore, the contributions of different 
sectors cancel each other, giving a vanishing density of topological charge at 8 = tt. 
In the thermodynamic limit, however, one sector could dominate (see below) and this 
would imply that the limits 8 — > tt and V — > oo do not commute, giving rise to the 
spontaneous Z 2 breaking. We cannot exclude that (x) Cjm vanishes, by accident, for 
some integer m at 8 = it. Usually however it is assumed that if the mean value of a 
local operator vanishes, it is because some symmetry forces it to vanish. In our case, 
as previously discussed, the "action" in the pseudo-partition functions Z c y(8 + 2irm) 
is not CP invariant at 8 = it. 

Equation ()2-7|) is valid for any value of 8. At 8 = (or at imaginary values of 8) 
and in the infinite volume limit the sector m = in (|2-7|) gives all the contribution 
to the partition function. This is a simple feature which follows from the fact that in 
this case all the terms in the summation (j2-4|) and the integrand in (|2-5[) are positive 
definite and then the partition functions (|2-4|) and (|2-5|) can be replaced, in the 
infinite volume limit, by the maximum of the summands and integrand respectively 
(saddle point solution). Since both maxima are coincident, we get the desired result. 

Even if we cannot exclude on theoretical grounds a phase transition in QCD at 
8 = 0, it seems not very likely. Then one expects that m = sector dominates at 
least in some interval around 8 = 0. If this sector dominates for every 8 between 
—it and 7r, we will get a phase transition at 8 = tt, with eventually a non- vanishing 
value of the topological charge density (remember the discussion following Eq. (|2-9|l ) 
i.e. the theory will show spontaneous CP breaking. If by accident the density of 
topological charge computed in the m = sector vanishes at 8 = tt, we should get 
a continuous phase transition since in any case sector m = — 1 will be the dominant 
one for 8 values between tt and 3tt. The only way to avoid a phase transition 
at 8 = tt would be that the vacuum energy density of the m = sector were a 
periodic function of 8 with period equal to 2tt. This of course is what happens 
in some trivial cases as QCD in the chiral limit, a limit in which there is no 8- 
dependence, or the dilute instanton gas approximation. 11 ' The last case is a model 
of non-interacting degrees of freedom, the partition function of which factorizes as the 
product of V identical partition functions; and it is well known that phase transitions 
and critical phenomena are collective phenomena which appear as a consequence of 
the interaction of infinite degrees of freedom. Systems with non-interacting degrees 
of freedom are equivalent to systems with one degree of freedom and therefore they 



9 -Vacuum 



5 



cannot show phase transitions. Therefore we can say that, excluding trivial cases, 
periodicity of the m = sector is not plausible. 

The other possibility is that m = sector dominates only until some critical 
9 C less than ir and then other sectors start to give a contribution to the partition 
function, in which case we will get a phase transition at this 9 C . 

Since in the previous argumentation we have not made use of any specific prop- 
erty of QCD, except the quantization of the topological charge, our result applies to 
any model with a quantized charge which appears as an imaginary contribution to 
the euclidean action. 



To see how this mechanism works in practical cases and to get intuition of 
what we can expect in physical systems, let us analyze in the following two simple 
examples: a model in which the probability distribution function of the density of 
topological charge at 9 = is assumed to be gaussian and the Ising model within an 
imaginary external magnetic field. 

The first example we want to discuss here includes models, as the quantum 
rotor with a density of topological charge at 9 = distributed according to a 
gaussian distribution. Thus, let us assume that the function hy(x n ) which enters 
Eq. (|2-4j) has the form 



where a is a parameter related to the width of the distribution. This form of hy(x n ) 
is a natural assumption from a physical point of view as a first approximation to 
the actual distribution of nearly any model. In fact, outside second order phase 
transitions, the probability distribution function of intensive operators as the density 
of topological charge is expected to be gaussian in the vicinity of its maximum. Of 
course, deviations from the gaussian behavior far from the maximum can induce 
important changes in the ^-dependence of the theory in the large 9 regime (as will 
become clear later). However, the gaussian distribution provides us with a simple 
model that can be analytically solved and gives useful insights on the general problem 
that we are addressing. 

The partition function of the model is 



Even if this model can be solved without the help of Eq. (|2-7|) . we will use it 
in order to clarify the relative weights of different sectors. The pseudo-partition 
function Z c y{9 + 2irm) entering Eq. (|2-7|) can be analytically computed, and the 
final result is: 



§3. Two simple examples 








(3-3) 
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Fig. 1. Expectation value of the topological charge density as a function of 8 in models with gaussian 
distribution. 

If \6\ < it, it is simple to verify that the free energy density f{9) is given by f(6) = 
limy^oo y log Zy{9) = —j^O 2 . The m = sector dominates for every 9 between 
— 7r and 7T. The vacuum expectation value of the density of topological charge, (x), 
is 



{X) =i Ya m<7T ' (3 ' 4) 

and the model breaks spontaneously parity at 9 = tt (see Fig. 1). 

This simple case that, as we have seen, illustrates very well the results obtained 
in the first part of the paper, has also further relevance. In fact, using a duality 
relation of gauge theories with a large number of colors and string theory on a 
certain space-time manifold, Witten has recently studied the 9 dependence of pure 
gauge theories in four dimensional* and found Eq. (|3-4|) , The combination of both 
results strongly suggests that the probability distribution function of the topological 
charge density in pure SU(N) gauge theory should be gaussian in the large N limit. 

The second illustrative example is the one-dimensional Ising model within an 
imaginary external magnetic field. The hamiltonian of this model can be written as 



H N = - J ^2 Si S i+ i - i — |— ^2 S i ' ( 3-5 ) 

i=l i=l 

where J is the coupling constant between nearest neighbors, kg is the Boltzmann 
constant, T the physical temperature, N the number of spins, and we assume periodic 
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boundary conditions. The partition function is given by 

Z N = ^^E^A+i+iOlE^ (3 . 6) 

{S} 

where F = J/ksT and the sum is over all spin configurations. 

For an even number of spins, the quantity 1/2 Y2i &i which appears in the imag- 
inary part of the hamiltonian is an integer taking values between —N/2 and N/2, 
and therefore it can be seen as a quantized charge. Furthermore, the theory has a 
Z2 symmetry at 9 = and 9 = it which, in the spirit of this work, is the analogue of 
CP in QCD. 

The transfer matrix technique allows to compute exactly the partition function 
defined in Eq. (pEK]i : Z N = A+ +A^, where A± are the two eigenvalues of the transfer 
matrix 

A±(0) = e F cos - ± (-e 2F sin 2 - + e~ 2F X * (3-7) 



2 V 2 
and we get for the mean value of the density of magnetization 

/ \ ■ sin(0/2) 



~ 4F - sin 2 



; (0/2)]V2' 



<vr. (3-8) 



The solution for every 9 between — tt and tt is the analytical continuation of the 
solution for a real magnetic field Since in the real magnetic field case the m = 
sector dominates always, this proves that the m = sector dominates for every 9 
between —tt and tt. A rather simple but long and straightforward calculation based 
on the analytical computation of the probability distribution function of the density 
of magnetization allows to demonstrate also this result, but the calculation is too 
long to be included here. 

For ^-values between tt and 3tt (—tt and — 3tt) the m = —1 (m = 1) sector 
dominates (A_ becomes larger than A+ in absolute value) and so on. Then we get, 
as expected, a periodic solution which shows a first order phase transition at 9 = tt 
with spontaneous Z2 breaking (see Fig. 2). 

In the ferromagnetic case (F > 0), the order parameter shows a divergence at a 
critical 9 C (T) running from to 71^'™ and the model is ill defined for 9 C (T) < 9 < tt. 



§4. Discussion 



Our main result, i.e., the existence of some singularity at 9 C between and 
tt, was conjectured for non-abelian gauge theories by 't HooflP in 1981. Based 
on the applicability of Poisson summation formula, it is verified by the two simple 
examples previously discussed and also by other models with well established results 
as CP N models in the strong coupling region,^ ^ two dimensional QED, 71 or 
high temperature QCD at imaginary chemical potentialp^^ where the imaginary 
chemical potential plays the role of 9 and the quantized charge is the baryonic charge. 
However it still seems to be in contradiction with some predictions of chiral effective 
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Fig. 2. Magnetization as a function of 9 in the one-dimensional antiferromagnetic Ising model. 



models of In these models the vacuum energy density is obtained, 

in the semi-classical approximation, as the solution of a given equation and the 
multiplicity of the solutions of this equation depends crucially on a certain inequality 
between quark masses. For realistic quark masses the solution would be unique, with 
periodicity 2tt, and the vacuum energy density would be an analytic function of 9 
for every 9, including 9 = 7r.^ 

There are however several possible explanations for the apparent contradiction 
between the chiral effective lagrangian results and ours, and we want just to enu- 
merate some of them here. First the chiral effective lagrangian does not belong to 
the general class of models to which our results should apply since the quantized 
charge does not appear in it as an imaginary contribution to the euclidean action. 
Second the vacuum energy density obtained from the chiral effective lagrangian is 
the classical solution, i.e., it does not take into account quantum fluctuations which 
could change the vacuum structure; and third, higher order corrections to the chiral 
effective lagrangian could also change the vacuum energy density even at the classical 
level (see for instance Ref.122))). 

We want to stress also that the results reported in this article are complementary 
to the general analysis on spontaneous CP breaking developed by Creutz in Ref. 23). 
In fact our results would exclude some general ansatzes for the effective potential 
compatible with the symmetries considered by Creutz in Ref. I23|) which drive to an 
order parameter analytical in 9. 
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